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Abstract 

We give sufficient conditions for the rigid body in the presence of an axisymmetric force field and 
a gyroscopic torque to admit a Hamilton-Poisson formulation. Even if by adding a gyroscopic torque 
we initially lose one of the conserved Casimirs, we recover another conservation law as a Casimir 
function for a modified Poisson structure. We apply this frame to several well known results in the 
literature. 

1 Introduction 

The generalized Euler-Poisson system, which is a dynamical system for the rotational motion of a rigid 
body in the presence of an axisymmetric force field, admits a Hamilton-Poisson formulation, where the 
bracket is the "— " Kirillov-Kostant-Souriau bracket on the dual Lie algebra e(3)*. The Hamiltonian 
function is of the type kinetic energy plus potential energy. The Hamiltonian function and the two 
Casimir functions for the K-K-S Poisson structure are conservation laws for the dynamic. 

Adding a certain type of gyroscopic torque the Hamiltonian and one of the Casimirs remain conserved 
along the solutions of the new dynamical system. In this paper we give a technique for finding a third 
conservation law in the presence of gyroscopic torque. 

We modify the K-K-S structure using the elements that generate the gyroscopic torque and we find 
necessary and sufficient conditions for this structure to be a Poisson structure. This leads us to what we 
call the Jacobi condition equation. Once this new Poisson structure is established, we give necessary 
and sufficient conditions for a function to be Casimir function. We call these conditions the Casimir 
condition equation. 

Using this framework, we study several general cases of gyroscopic torques, that include some well 
known situations found in the literature. We encompass the gyrostatic torques, affine gyroscopic torques, 
nonlinear torques studied by Yehia pQ , [5] , gyroscopic torques along one axis of inertia studied by Borisov 
and Mamaev [3J. 

It would be interesting to relate the methods presented in this paper with the methods of constructing 
and analyzing gyroscopic torques given in Chang and Marsden [4] and in Chang et al. [5]. The 
understanding of the link between the two methods could give better results on studying the dynamics 
in the presence of gyroscopic torque. 

2 Mechanical aspects 

The equation for the rotation of a rigid body about a point O is 

luj + Q x Ilo = M 
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in the cases: 

(i) the point O is fixed with respect to an inertial frame; 

(ii) the point O coincides with the center of mass. 

In the above equation lu is the angular velocity, I is the moment of inertia tensor at the point O and 
Mo is the applied torque about the point O. For the following considerations we use the body frame H 
with the origin O and its axes are principal axes of inertia of the rigid body at the point O. The matrix 
of the moment of inertia tensor in this body frame has the form 

fa =diag(A,B,C), 

where A, B and C are principal moments of inertia. 

Our rigid body is situated in an axisymmetric force field. We denote by 7 a vector which generates 
the axis of symmetry and we have the Poisson's differential equation 

7 + LU x 7 = 0, 

with the time derivative taken with respect to the body frame 3?. 

The rigid body is acted, also, by a gyroscopic torque. We study the case in which the applied torque 
Mo has the form 

M = M a + M 9 , 

where M a is the torque generated by the axisymmetric force field and M ff is the gyroscopic torque. The 
torque M a depends on 7 and in a sufiiciently extensive class of cases it is represented in the form 

M a = 7 x V 7 Z7, 

where U is a smooth function which depends on 7. The most general form of the gyroscopic torque 
considered in this paper is 



M„ 



-LU X /I, 



(2.1) 



where pi is a vectorial function which depends on Co and 7. 

For mathematical reason we use the angular momentum vector M to describe our problem. Its 
connection with the angular velocity is given by 

M = leu. 

According with the above considerations the system which describes the rotational motion of a rigid 
body situated in an axisymmetric force field and in the presence of a gyroscopic torque has the form 



M = -r'M x (M + f!(M, 7)) + 7 x VyUtf) 
1 



7 = 7 x Ji _1 M. 

In the absence of the gyroscopic torque, i.e. /2 = 0, the system 



(2.2) 



becomes 



M = -IT 1 !/ x M + 7 x V 7 C/(7) 
7* = 7 x l- x M 



(2.3) 



and we obtain what is known as the generalized Euler-Poisson system. For /2 = and U = mg^ ■ 7 one 
obtains the classical Eulcr-Poisson problem which describes (on the leaf ||7|| = 1) the rotational motion 
of the rigid body under gravity. The constant m is the mass of the rigid body, g is the acceleration of 
gravity and £ is the vector of position of the mass center. 

This system has the following Hamilton-Poisson structure. The Poisson bracket is the " — " Kirillov- 
Kostant-Souriau bracket on the dual Lie algebra e(3)*. In the coordinates (Mi, M%, M3, 71, 72, 73) it 
has the expression 
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The Casimir functions are 

C 1 (M,^) = ^f 

and 

C 2 (M,7)=M-7. 

The Hamiltonian function is 

H(M,i) = -M-l- 1 M + U(j). (2.4) 

In the presence of a gyroscopic torque of the form (|2.1[) or equivalently M g = — l~ l M x jl the above 
functions H and C\ remain constants of motion along the solutions of (|2.2I) . The time derivative of C2 
along the solutions of (12.21) is 

^C 2 =7-(/IxI- 1 M). 

Unless /I is a linear combination of 7 and E M = u, the function C2 is not a constant of motion for 
the system (ET2"j) . 



3 Hamilton-Poisson formulation in the presence of gyroscopic 
torque 

As we have seen in the previous section, when we add a gyroscopic torque to the Hamilton-Poisson 
system (|2.3[) we lose one of the Casimirs as a constant of motion for the system (|2.2p . 

We will recover a third constant of motion as a Casimir function by modifying the K-K-S Poisson 
structure. Our modification of the Poisson structure will be such that the system (|2.2j) becomes a 
Hamilton-Poisson system. We consider the following skew-symmetric 6x6 matrix 
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where the components ^i,/i2,/^3 of the vector field jl are smooth functions of the variables M and 7. 
The matrix II^ generates a Poisson structure if and only if the Jacobi identity is verified, i.e. 



0. 



for all i,j,k = 1,6. The above system of partial differential equations reduces, in our case, to the 
following system of ten PDEs for the functions /ii , ^2 , M3 : 



7XV A jMfc = 0, fc = l,2,3 

7 ■ curLy/2 + M ■ curl A |/i + /1 ■ curl A i/i 



0. 



(3.1) 



Solving independently the systems of equations 7 x V^j/ifc = for k = 1, 2, 3 we obtain /Zfe as 
functions of four variables 71, 72, 73 and s := M ■ 7, 

Mfe = Vk{l,M ■ 7). 



Substituting in the last equation of the system (|3.ip the /ifcS found above we obtain the Jacobi condi- 
tion equation for //: 

7 • curLy/2 + ^ ■ (7 x 94/2) = 0. (3.2) 

The partial derivative 84 jl is the derivative with respect to the variable s = M ■ 7. If jl verifies (|3.2|) . 
then 11^ generates a Poisson bracket on R 6 . 
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A smooth function C(M,j) is a Casimir function for the Poisson structure II^ if and only if it 
satisfies the following system of six partial differential equations: 



f 7 x V a C = 

\ (M + ft) x V a5 C* + 7 x V^C* = 0. V ' ' 

Solving the system 7 x V A ^C = we obtain C as a function of four variables 71,72,73, s = M ■ 7, 
i.e. C = C(fy, M • 7). Substituting C found above in the system (M + ft) X V^C + 7 x V^C = we 
obtain the Casimir condition equation: 

7 x ((d 4 C)fl - V^C*) = 0. (3.4) 

We note that C\(M, 7) = -7 2 remains a Casimir for 11^ as it verifies the Casimir condition equation 
(|3.4j) . Next, we summarize the considerations we have made in this section. 

Theorem 3.1. The following results hold: 

(i) Up generates a Poisson bracket {•, -}p if and only if the vectorial function fx depends only on the 
four variables 71,72, 73, M • 7 and it verifies the Jacobi condition equation (|3.2[) . 

(ii) A smooth function C(M, 7) is a Casimir function for the Poisson bracket {•, -}p if and only if C 
is a function that depends only on the four variables 71,72,73, M ■ 7 and it verifies the Casimir 
condition equation (|3.4|) . 

(Hi) If Up generates a Poisson bracket, then the system (|2.2j) has the Hamilton- Poisson formulation 
(R 6 , {•, -}p, H), where H is the Hamiltonian function (|2.4[) . 

For a gyroscopic torque that verifies the conditions of the above theorem we provide an algorithm 
that allows us to recover a third constant of motion that we have lost when we added the gyroscopic 
torque to the generalized Euler- Poisson system (|2.3|) . This third constant of motion will be a Casimir 
function for {•, -}p, independent of the Casimir G\. 

4 Applications 

We will discuss several well known gyroscopic torques that appear in the literature. We will study these 
examples in the theoretical frame that we have developed in the previous section. 

4.1 The case jl = jj(^f) 

In this particular case when jl depends only on the Poisson variables 71,72,73 the Jacobi condition 
equation (|3.2p becomes 

7 • curl^/I = 0, (4.1) 
which is equivalent with the vector field 7 x /2(7*) being a solenoidal field, i.e. 

div^(7 x ft) = 0. 

A large class of examples corresponds to the following solution for equation (14. ip . 

Al) =^(7)7 + V^(7), (4.2) 

where tp,tp are smooth functions depending only on variables 71,72,73- 
For the above jl the Casimir condition equation (|3.4[) becomes 

7 x ((d 4 C)(V(707 + V^(7*)) - V^C) = 0. 

We search a solution under the hypothesis d^C = 1. Consequently, we have the following Casimir 
condition equation, 

7 x V=f(ip - C) = 0. 
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A particular solution is given by 

C(7,M-7) =M -7+^(7). (4.3) 

1 -,2 

This Casimir is functionally independent of the Casimir C\ = — 7 . 
In Yehia [I] it is considered a torque of the form 

Al) = -(div/)7 + V^(f- 7), 

where / = £(7) is a smooth function. Choosing ^(7) = — div^Z and </?(7) = / -7, jl can be put in the form 
(|4.2|) and consequently these torque cases admit a Hamilton-Poisson formulation. The second Casimir 
is given by C(M,"f) = (M + 1) • 7. This Casimir appears in [T] as the cyclic integral of motion. 

The case of gyrostatic torque 

An important case is the gyrostatic torque M g , where ^(7) = ^0 is a constant vector. A gyrostatic 
torque can be produced by mechanical rotors or it is the consequence of the liquid action which is 
contained in rigid body cavities. An interesting situation for which we have a gyrostatic torque is 
presented in the Zhukovsky theorem (see [B], pp. 51), where a fluid mass with an initial velocity in a 
multiply connected cavity also performs an action that is similar to the action of some rotor attached 
to the rigid body (the fluid completely fills the cavities of the rigid body and the motion of the fluid is 
potential). Various mathematical aspects, for particular functions U, are studied in some papers. For a 
linear function U, Gavrilov investigates, in [7], the integrability of the system. 

A gyrostatic torque can be put in the form (|4.2I) by choosing ^(7) =0 and <p(j) — flo-J- From (|4.3|) 
we obtain the second Casimir 

C(M,1) = (M + p )-r ▲ 

The case of affine gyroscopic torque 

Grioli studies in [8] the case of the rotation of a charged rigid body, with a fixed point, situated in 
a constant force field. We present some details of the calculation of the torque generated by this force 
field. We denote by O the fixed point, q the electric charge, 7 the magnetic field and D the domain of 
the rigid body. Let P a particle of the rigid body and v(P) its velocity. The particle is acted by the 
Lorenz force 

F = q(v x 7) = q((uj x OP) x 7). 
The torque generated by the magnetic field is 

Mo = / q(P)OP x Fdv = uj x £7, 
Jd 

where the tensor £ is defined by 

£ = / q{P)OP <g> OPdv. 
Jd 

It is easy to see that the torque generated by the magnetic field is a gyroscopic torque of the form (|2.ip 
with [1 = £7. The matrix associated with the tensor £ is a 3 x 3 constant symmetric matrix with respect 
to the body frame 51. In [5] it is considered a generalized Grioli problem for which the gyroscopic torque 
is generated by fl = £7 + e, where e is a constant vector in the body frame 3?. 

In our case we consider jKfj) — A7 + /lo, where A is a 3 x 3 constant matrix and /2o is a constant 
vector, both computed in the body frame. It is easy to verify that such a torque /2(7*) satisfies the 
Jacobi condition equation if and only if the matrix A is symmetric. In this case jl can be but in the 

form (|4.2|) by choosing ^(7) = and ^(7) = -7 • Aj + ft ■ 7. 

Wc obtain the second Casimir C(M,j) — -7 ■ /I7 + (M + jlo) ■ 7. Consequently, the generalized 
Grioli problem admits a Hamilton-Poisson formulation. ▲ 

Nonlinear cases 

In what follows we will analyze several particular cases of systems of type (|2.2[) which are presented 
in Yehia [I] and [2]. 
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a) Wc consider a rigid body with the Kovalevskaya configuration (A = B = 2C), where the poten- 
tial is given by 

U = C(ai7i + a 2 j 2 ) - Ckj 3 (n + rii 7 i + n 2 -f 2 ) - \ C ( n + n i7i + ™272) 2 (27 2 + 272 + I3) 

and the gyroscopic torque is generated by 

/ii = C(-nji - ni7i + 2 «i72 + n ili ~ 3«27i72) 
M2 = C{-ny 2 + 2n 2 7i - n 2 j 2 + n 27 f - 3ni7i7 2 ) 
jU 3 = C(k - 37173 - 5ni7i7 3 - 5n 2 7273), 

where 01, a 2 , k, n, rii, n 2 are real constants. 

This gyroscopic torque can be put in the form (|4.2[) with 

tp = C[fc 73 + (n + ni7i + n 2 72)(2 7l 2 + 2 7 | + 7 f )] 

and 

^ = — C(5n + 7rti7i + 7n 2 j 2 ). 
Consequently, this case has a Hamilton-Poisson formulation and the second Casimir is given by 

C(M, 7) = M • 7 + C[fc 73 + (n + 71x71 + n 2 y 2 )(2jl + 2 7 2 + 7 f )], 

which coincides with the cyclic integral of motion _?2 from [T], pp. 342, subsection 4.1.1. <> 

b) We consider a rigid body with the Kovalevskaya configuration (A = B = 2C), where the poten- 
tial is given by 

U = C(ai7i + 0272) + , = 5 - \c[n+ n l7l + n 2l2 + = I (2 7l 2 + 2 7 | + 73 2 ) 
V7i+7 2 2 V V7f+75/ 

and the gyroscopic torque is generated by 

Mi = G ( -n 7i - ni7 2 + 2rai 7 f + ni7 2 - 3n 2 7i7 2 + ' ' ^ 



(7i 2 + 7 2 2 )^ 

M2 = C [ -n 72 + 2n 2 7 2 - n 2 72 + "273 ~ 3«i7i72 + . J^ 12 ^ 3 
V ~ (7f + 7 2 ) 2 

( Nj 3 \ 

/j 3 = -C73 3n + 5ni 7 i + 5n 2 72 + 



VtF+tI/ 

where 01, 02, £, N, n, rii, n 2 are real constants. 

In Yehia [I], pp. 343, subsection 4.1.2, the function 

h = M-j + C[n + n l7l + n 27 2 + , - , I (2 7 i 2 + 2 72 2 + 7a) 
\ V7f+72/ 

is asserted to be an integral of motion. By direct computation it can be verified that this statement 
does not hold. We will present a method to solve this incompatibility. 

If we want to keep the gyroscopic torque as above, then we observe that it verifies the Jacobi 
condition equation (|4.1|) and from the Casimir condition equation (|3.4j) we obtain 



C(M, 7) = M • 7 + C (n + n l7l -I- 71272) (2 7l 2 + 2 7 | + 7 |)+ 

+ l ( 73(7^71) -2 _ (7? + 7| + 7|) arctan . 73 



2 v v^r^f x/Tf+Tf 



G 



If we want to keep J 2 as a constant of motion we will find a new gyroscopic torque generated by 
p' = (pi, fj,2, M3), which can be put in the form (|4.2p and for which 7 2 becomes a Casimir function for 
the Poisson bracket generated by p' . We choose 

<P(Tf) = c(n + mji + n 272 + —?L^ J (2 7l 2 + 2 7 | + 73 2 ) 

and 

4>(l) = -C\hn + 7ni71 + 7n 272 + TV ■ 2(71 + J 2) ~ 7 3 3 ~ 1 ) . 

V (7i+7a) 2 / 

We observe that /xi = ^(7)71 + ^§^-, = ^(7)72 + an d consequently /i 3 is given by 

, gg(j) 7 3 (7i 2 +72 2 )+7 3 2 + l 
M3 = W 73 + = M3 + CW73 — —3 

a 73 ( 7 i 2 + 7 2 2 )^ 

= -C 73 ( 3n + SniTi + 5n 272 - AT - 7 3 + 9 \ 3 ) ■ 
V (7i 2 +7l)'/ 

The first, third and fourth cases from Yehia [2J can be analyzed as in the case a); it is also the 
situation for the generalization of Lagrange case from Yehia [2J. The third case from Yehia [I], the 
second and fifth cases from Yehia [2J can be analyzed as in the case b) . ▲ 

4.2 The case ft = #(7, M ■ 7) 

In this section we will study two general subcases of gyroscopic torques that depend also on s — M ■ 7. 
We will find a Hamilton-Poisson formulation for these subcases. 

The subcase p = a(s)V^</?( 7 ) + 6(7, s)j 

For this subcase we consider a:/Cl4l, 6:Dx/CR 3 xt^land(^:LiCl 3 ^K smooth 
functions, where I and D are domains in the corresponding spaces. 

The gyroscopic torque generated by p verifies the Jacobi condition equation. Indeed, 7 -curl^/i = 
and 



P- (7 X dip) = -7* • (p X 9 4 m) = 7 ' ( a(s)^(7>)V^( 7 ) x 7 + 6( 7 ,s)— 7 x V^( 7 ) ) = 0. 



The Casimir condition equation becomes 

7 x ((9 4 C)a(s)V^( 7 ) - V^C) = 0. 

We will search for Casimir functions of the form 

Ctf,s) = D(j)E(s). 

Introducing this expression for the Casimir in the above equation we obtain, 

7 x (E'( S )a(s)D(i)\7 W (i) - E(s)^D(j)) = 0. 

Sufficient conditions for the above equation to be satisfied are that the functions E and D verify the 
following equations: 

a(s)E'(s) = E(s) 

and 

£>( 7 )V^( 7 ) = V ? D( 7 ). 
Consequently, we obtain a particular solution 



and 

D{l) = exp^O?). 
Thus, a Casimir function for 11^ is given by 

C(1,s)= [ -L-ds + tptf). (4.4) 
J a{s) 

This case is a generalization of the gyroscopic torque given by relation l|4.2[) for a(s) = 1 and 
6(7, s) = The Casimir (I4.3[) is a particular case of the Casimir function (|4.4I) . ▲ 

Gyroscopic torques along one axis of inertia 

Without losing the generality we will study gyroscopic torques generated by /2 = (0, 0, ^3(7, s)). 
The Jacobi condition equation for /2 of this form reduces to 

71 o 72^— =0, (4.5 

072 071 

which has the solution /x 3 = ^(r, 73, s), where r := ^(Ti + 7a )• 

The Casimir condition equation is equivalent with the following system of equations: 

„, /ac„ acA , „, ac _ n 
72 (-57^3 - ^ ) + 73 - 

ry, ( Q£ _ ac\ , ac = 

From the first equation we obtain that the Casimir function has the form 

C = C(r, 73 ,s). 

Substituting this form in the last two equations we obtain 

dC dC dC n fl „. 
^TT -^7-+737r=0. (4.6) 
us c/73 or 

For /X3 of the particular form 

^3(^,73,5) = /8(7 3 )<J(s) 

we have a Casimir function given by 

C(r,73,s) = J ^ds + J f3(j 3 )dj 3 



as a solution of 

In Borisov and Mamaev [3] it is introduced a gyroscopic torque generated by H\ = [12 = and 
M • 7 

/13 = for the Kovalevskaya configuration of the rigid body and the potential U — aOy? —72)1 

73 

where a is a constant. This choice for fl verifies (|4.5[) and the function C(r, 73, M ■ 7) = 73s = 73 (M ■ 7) 
verifies (|4.6|) and thus is a Casimir function for Hp. ▲ 
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